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Abstract 

We study the electronic energy spectra and wave functions on the square Fibonacci 
tiling, using an off-diagonal tight-binding model, in order to determine the exact nature 
of the transitions between different spectral behaviors, as well as the scaling of the total 
bandwidth as it becomes finite. The macroscopic degeneracy of certain energy values in 
the spectrum is invoked as a possible mechanism for the emergence of extended electronic 
Bloch wave functions as the dimension changes from one to two. 

1 Background and Motivation 

We continue our initial studies [1] of the off-diagonal tight-binding hamiltonian on the square 
Fibonacci tiling [2], in order to gain a better quantitative understanding of the nature of the 
transitions between different spectral behaviors in this 2-dimensional (2d) quasicrystal. We also 
consider more carefully the transition of the spectrum from singular-continuous to absolutely 
continuous, in going from one to two dimensions, and the implications of this transition on the 
possible emergence of extended Bloch wave functions. 

The square Fibonacci tiling is constructed by taking two identical grids — each consisting 
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of an infinite set of lines whose inter-line spacings follow the well-known Fibonacci sequence of 
short (S) and long (L) distances — and superimposing them at a 90° angle. This construction 
can be generalized, of course, to any quasiperiodic sequence as well as to higher dimensions. 
On this tiling we define an off-diagonal tight-binding hamiltonian with zero on-site energy, and 
hopping amplitudes: t for vertices connected by short (S) edges; 1 for vertices connected by 
long (L) edges; and zero for vertices that are not connected by edges. 
The 2-dimensional Schrodinger equation for this model is given by 

t n+1 *£>(n + 1, m) + t n V(n - 1, m) 

+t m+1 ^(n, m + 1) + i m *(n, m - 1) = EV(n, m), (1) 

where ty(n, m) is the value of a 2d eigenfunction on a vertex labeled by the two integers n 
and to, and E is the corresponding eigenvalue. The hopping amplitudes tj are equal to 1 or t 
according to the Fibonacci sequence as described above, and we do not lose much by limiting 
ourselves to the case t > 1, as the important features of the model depend on the ratio of the 
two hopping amplitudes. 

With no additional assumptions other than the absence of diagonal hopping this 2-dimensional 
eigenvalue problem, as well as its higher dimensional versions, are all separable. This allows us 
to use the known solutions for the Id problem [3, 4] in order to construct the solutions in two 
and higher dimensions (as was done for similar models in the past [5, 6]). Two-dimensional 
eigenfunctions can therefore be expressed as products 

^ij(n, m) = t/ji(n)ipj(m), (2) 

with eigenvalues 

Eij = Ei + Ej, (3) 

where ipi(n) and ipj(n) are two of the eigenfunctions of the corresponding Id eigenvalue equation 
on the Fibonacci chain, with eigenvalues Ei and Ej. Consequently, the 2d density of states 
(DOS) is given by a convolution of two Id densities of states. 
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2 Study of the 2d spectrum 



Although the spectrum of the Id model is always a Cantor-like set with zero bandwidth and 
an infinite number of bands, the 2d spectrum exhibits different behavior for different ranges 
of the parameter t. For t close to 1 the gaps in the Id spectrum are small and disappear as 
a result of the addition of spectra (3), producing an absolutely continuous 2d spectrum. For 
i > 1 the 2d spectrum retains the Cantor-like structure of the Id spectrum. For intermediate 
values of t the number of bands tends to infinity but the total bandwidth remains finite and it 
is unclear whether the spectrum has an absolutely continuous part. 

In order to study the exact structure of the 2d spectrum — the total bandwidth and number 
of bands as a function of t — one has to perform a convolution of two Id spectra. The numerical 
study of the Id spectrum is difficult for low values of t (t s=s 1) for which the gaps are very 
small, and for t 3> 1 for which the bands are very small. Further numerical difficulties are 
encountered in the convolution process itself. The desired 2d spectrum, consisting of the set of 
points for which the 2d DOS is not zero, is given by 



and Fn is the N th Fibonacci number. Because Fn ~ t n the number of intervals in the 
union (4) behaves like t 2N , and although many of the intervals may overlap, any attempt to 
directly study the band structure and the DOS requires a construction of the whole set. The 
amount of memory required to store such sets makes exact calculation implausible for high 
values of N regardless of the value of t. Luckily, the total bandwidth of the 2d spectrum can be 
calculated without full knowledge of the DOS. This is done by examination of a set of uniformly 
distributed energies in the relevant range and testing for each energy whether it lies within one 
of the above intervals. For a sufficiently large set of energies this procedure yields the measure 
of the spectrum after proper normalization. 

For given values of N and t the Id spectrum is a symmetric set of Fn intervals within 




(4) 



where a, 



is the i th interval in the relevant Id spectrum of the N th approximant, 
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Figure 1: Total bandwidth of the 2d spectrum divided by 4rj^(t). For values of t up to ~ 2 
the spectrum consists of a single band regardless of N. At higher values the bandwidth clearly 
depends on N. 



[— r N(t), rjv(t)], where r^(t) is the maximal eigenenergy of the N Id approximant. The 2d 
spectrum is therefore a subset of the interval [— 2r^{t), 2r7v(t)], which gives an upper bound of 
^ r N{t) on its total bandwidth. Fig. 1 shows the total bandwidth divided by 4r7\r(£) as a function 
of t for a sequence of approximants iV = 6 ... 17. For values of t up to ~ 2 the spectrum 
consists of a single band. For t's above 2 the normalized bandwidth becomes dependent on 
N. Using finite-size scaling (as suggested by Ashraff et al. [6]) we can obtain an estimate for 
the critical value t c at the transition between the intermediate and large t regimes, where the 
total bandwidth vanishes. Fig. 2 shows the total bandwidth as a function of 1/N for different 
values of t > 2. Though some oscillatory behavior is observed, the behavior seems quite linear 
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Figure 2: Total bandwidth as function of 1/N for different values of t. Linear extrapolation to 
1/N — > gives a prediction for the total bandwidth at the infinite N limit. 

and extrapolation to 1/N — » yields a prediction for the total bandwidth in the infinite iV 
limit. By plotting these estimates as a function of t we obtain Fig. 3 which again exhibits an 
overall linear behavior. Extrapolation yields an estimate for a critical value of t c ~ 5.375 for 
the transition between the intermediate and large t regimes. 

It is worth noting that as t is increased and the gaps grow, the degeneracy of energy levels 
decreases. At the asymptotic limit each interval in the 2d spectrum should be exactly doubly 
degenerate (neglecting a single interval around E = 0). When this limit is obtained the 2d 
bandwidth should be approximately equal to F^B^t, N), where Bid(t,N) is the bandwidth 
of the corresponding Id spectrum. The numerical data shows that this limit is to be obtained 
well beyond t c . It is hence evident that the large t regime of the 2d spectrum is further divided 
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Figure 3: Extrapolation of the predicted infinite limit bandwidth to find the critical value t c , 
above which the total bandwidth tends to zero. 

into two regimes, one (at extremely large t) in which the energy levels are no-more than doubly 
degenerate, and another (at moderately large t) in which energy levels may have a higher 
degeneracy. 

3 Emergence of extended wave functions 

Since the eigenfunctions of the Id hamiltonian are critically localized it is clear that the 2d 
eigenfunctions (2) should also be critically localized (Fig. 4(a-c)). However, the existence 
of a regime for t close to 1 in which the spectrum, or at least a part of it, is absolutely 
continuous hints at the existence of extended wave functions. We proposed in the past [1] that 
the degeneracy of the 2d spectrum, at values of t close to 1, may provide an explanation for 
the emergence of extended wave functions as the dimensionality increases. 

To illustrate the possible role of degeneracy in the emergence of extended wave functions 
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Figure 4: Wave functions on a single 144 x 144 2d approximant. (a) A critical wave function 
^n,m', (b) A symmetrized critical wave function ^= (^ n ,m + *m,n); ( c ) A zero energy wave 
function \&_ njn ; and (d) A linear combination of zero energy wave functions producing an 
extended Bloch wave function along the diagonal, clearly showing the underlying Fibonacci 
modulation. 
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we analyze the E = eigenfunction. These functions are macroscopically-degenerate for any 
value of t. A set of Id eigenfunctions can be found by directly diagonalizing the Id hamiltonian 
for a single approximant with periodic boundary conditions. The F/v eigenfunctions obtained 
in this way generate eigenfunctions for the 2d model, Fn of which have E = 0. These 
eigenfunctions, given by *$>i(n,m) = (—l) m "il>i(n)'ipi(m), are each critically localized, but we 
can demand that a linear combination YliLi a i^i{ n i m ) °f them have certain values on a given 
subset of sites by solving the equations for the coefficients a^. Setting the values of a linear 
combination to be identical on half the points along the diagonal results in eigenfunctions of 
the form shown in Fig. 4(d). This wave-function is an extended Bloch function along the 
diagonal — clearly showing the underlying Fibonacci modulation — but is strongly localized in 
the perpendicular direction. 

Whether the degeneracy of other energy values is sufficiently high to allow for a similar 
process to be carried out is not clear, but the chances for success should strongly improve as 
the tunnelling amplitude t approaches 1. This also corresponds to the regime in which the 
spectrum is (at least partly) absolutely continuous. For large t, where the spectrum is singular 
continuous, the degeneracy of levels other than E = is very small and no extended wave 
functions are expected to be found. Note that in 3 dimensions all the energies of the original Id 
spectrum are macroscopically degenerate. This means that extended wave-functions, similar to 
that of Fig. 4(d), can be constructed for a wider range of energies. Thus, as the dimensionality 
further increases so does the degeneracy of energy levels, and with it the ability to construct 
Bloch wave functions. 
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